vAOMList = [4.6, 3.2, 2.0, 1.7, 1.4, 1.0, 0.7];
powerList = [150, 100, 45, 30, 20, 10, 5];
initBx0List = [-0.01, 0.005,  0.02, 0.025,  0.03,  0.03,   0.04];
initBy0List = [-1.28, -1.30, -1.32, -1.32, -1.32, -1.32, -1.325];

for k = 1:length(vAOMList)
    nmrg.GroupCharacterization( 110,  powerList(k),  1.8,  377.10280, 384.21000, ...
        'initialField', 5.9300, 'isFixBz0', true, 'Xe129Coeff', 253.5/5.93, ...
        'AOMVoltage', vAOMList(k), 'EaterVoltage', 0.4, 'initialBx0', initBx0List(k), 'initialBy0', initBy0List(k), 'amplitudeTransX', 0.01, 'amplitudeTransY', 0.05, ...
        'Xe129DrivingAmplitude', 0.016, 'Xe131DrivingAmplitude', 0.050, 'EstXe129T2', 20.0, 'EstXe131T2', 30.0, 'EstXe129T1', 180.0, 'EstXe131T1', 60.0, ...
        'description', sprintf('Power%dmW', powerList(k)));
    pause(10);
end